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Abstract. In this work we study the internal spatial structure of 16 open clusters in the Milky 
Way spanning a wide range of ages. For this, we use the minimum spanning tree method (the 
Q parameter, which enables one to classify the star distribution as either radially or fractally 
clustered), King profile fitting, and the correlation dimension (D c ) for those clusters with fractal 
patterns. On average, clusters with fractal-like structure are younger than those exhibiting radial 
star density profiles. There is a significant correlation between Q and the cluster age measured 
in crossing time units. For fractal clusters there is a significant correlation between the fractal 
dimension and age. These results support the idea that stars in new-born clusters likely follow 
the fractal patterns of their parent molecular clouds, and eventually evolve toward more centrally 
concentrated structures. However, there can exist stellar clusters as old as ~ 100 Myr that have 
not totally destroyed their fractal structure. Finally, we have found the intriguing result that the 
lowest fractal dimensions obtained for the open clusters seem to be considerably smaller than 
the average value measured in galactic molecular cloud complexes. 
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1. Introduction 

The hierarchical structure observed in some open clusters is presumably a consequence 
of its formation in a turbulent medium with an underlying fractal structure (Elmegreen 
& Scalo 2004). Otherwise, open clusters having central star concentrations with radial 
star density profiles likely reflect the dominant role of gravity, either on the primordial 
gas structure or as a result of a rapid evolution from a more structured state (Lada & 
Lada 2003). Therefore, the analysis of the distribution of stars may yield information 
on the formation process and early evolution of open clusters. It is necessary, however, 
that this kind of analysis is done by measuring the cluster structure in an objective, 
quantitative, as well as systematic way. Here we study the internal spatial structure in a 
sample of 16 open clusters spanning a wide range of ages. 

2. Procedure 

(a) We first used VizieR (Ochsenbein et al. 2000) to search for catalogs containing 
both positions and proper motions of stars in open cluster regions. 

(b) We applied a robust non-parametric method to assign cluster memberships (Cabrera- 
Caho & Alfaro 1990). This method makes no a priori assumptions about cluster and field 
star distributions. 

(c) We fitted King (1962) profiles to the radial density distribution of cluster members. 
From these fits we obtained both the core radius (R c ) and the tidal radius (i?t). 
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Figure 1. The minimum spanning tree is the set of straight lines connecting the points such 
that the sum of their lengths is a minimum. Here we show minimum spanning trees for three 
open clusters, from which we can calculate the structure parameter Q. Star positions are indi- 
cated with blue circles and red lines represent the tree. The value of Q quantifies the spatial 
distribution of stars. For IC 2391 the stars are distributed following an irregular fractal pattern 
(Q = 0.77 < 0.8), for M 34 the stars are distributed roughly homogeneously (Q = 0.8), and for 
M 11 the stars follow a radial density profile (Q — 1.02 > 0.8). 

(d) Then, we used the minimum spanning tree technique (see Fig.QJ to calculate the 
dimensionless parameter Q (see details in Cartwright & Whitworth 2004 and Schmeja & 
Klessen 2006). The value Q ~ 0.8 separates radial clustering (Q > 0.8) from fractal type 
clustering (Q < 0.8). 

(e) Finally, we calculated the correlation dimension (D c ) and its associated uncertainty 
by using an algorithm which gives reliable results (Sanchez et al. 2007a, Sanchez & Alfaro 
2008). 

3. Main results 

Table [T] summarizes the relevant data (ages and distances were taken from the Webda 
database). On average, stars in young clusters tend to be distributed following clustered, 
fractal- like patterns (Q < 0.8), whereas older clusters tend to exhibit radial star density 
profiles (Q > 0.8). However, the statistical analysis indicates that there is no significant 
correlation between Q and log(T). If instead we consider the variable T/R t , which is 
proportional to the cluster age measured in crossing time units (assuming nearly the 
same typical velocity dispersion for the open clusters), then a significant correlation 
is observed (Fig. [5]). Additionally, we observe significant correlations (confidence levels 
above 96 %) between D c and T (cluster age) and also between D c and T/R t (age in 
crossing time units) for those clusters with internal substructure (Fig. [3]). 

4. Discussion 

Our results support the idea that stars in new-born cluster likely follow the fractal 
patterns of their parent molecular clouds, and that eventually evolve toward more cen- 
trally concentrated structures (see Schemja & Klessen 2006; Schmeja et al. 2008, 2009; 
Sanchez et al. 2007a, 2009). However, this seems to be only an overall trend. The very 
young cluster a Orionis (age ~ 3 Myr) exhibits a radial density gradient with Q ~ 0.88 
(Caballero 2008). On the other hand, Table [1] shows open clusters as old as ~ 100 Myr 
that have not totally destroyed their clumpy structure (for example, both NGC 1513 and 
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Table 1. Properties of the clusters 


in the sample. 






Name 


log T 
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D c 


IC 2391 
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1 1 


DZ 


1 AG. 
1.40 


2.65 
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1.74 ±0.20 


M 11 
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NGC 188 




OA/I 7 

ZU4 ( 


1 /I CCA 

i4oy 


O AA 
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0.91 




NGC 581 


1 .336 
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o2o 


1.38 


11.86 


0.76 


1.79 ±0.06 


NGC 1513 
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1 QOA 


lOO 


1.00 


7.73 


0.72 


1.82 ±0.09 


NGC 1647 
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1.94 ±0.02 


NGC 1817 
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Z II 
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11.97 


0.79 


1.94 ±0.04 


NGC 1960 


7.468 


1318 


311 


2.96 


8.77 


0.87 




NGC 2194 


8.515 


3781 


228 


3.17 


10.31 


0.85 




NGC 2548 


8.557 


769 


168 


2.61 


9.16 
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NGC 4103 


7.393 
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0.72 
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1.85 ±0.04 
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NGC 5281 


7.146 


1108 


80 


0.62 


2.44 


0.84 




NGC 6530 


6.867 


1330 


145 


1.43 


7.47 


0.67 


1.74 ±0.09 



Note: T = cluster age (Myr); D = distance (pc); N s = number of members; R c = core radius 
(pc); Rt = tidal radius (pc); Q = structure parameter; D c — correlation dimension. 



NGC 1647 have Q ~ 0.7). Goodwin & Whitworth (2004) simulated the dynamical evo- 
lution of young clusters and showed that the survival of the initial substructure depends 
strongly on the initial velocity dispersion. Fractal clusters with a low velocity dispersion 
tend to erase their substructure rather quickly. However, if the velocity dispersion is high, 
such that the cluster remains supported against its own gravity or even expands, then 
significant levels of substructure can survive for several crossing times. Thus, our results 
give some observational support to Goodwin & Whitworth's (2004) simulations. 

From Fig. [3l we can see that clusters with the smallest correlation dimensions (D c = 
1.74) would have three-dimensional fractal dimensions around Df ~ 2.0 (estimated from 
previous papers, see e.g. Fig. 1 in Sanchez & Alfaro 2008). This is a very interesting result 
because this value is considerably smaller than the average value estimated for galactic 
molecular clouds in recent studies, which is Df ~ 2.6 — 2.7 (Sanchez et al. 2005, Sanchez 
et al. 2007b). Young, new-born stars probably will reflect the conditions of the inter- 
stellar medium from which they were formed. Therefore, a group of stars born from the 
same cloud, i.e. born at almost the same place and time, should have a fractal dimension 
similar to that of the parent cloud. If the fractal dimension of the interstellar medium has 
a nearly universal value around 2.6-2.7, then how can some clusters exhibit such small 
fractal dimension values? Perhaps some clusters may develop some kind of substructure 
starting from an initially more homogeneous state. This possibility has been confirmed 
in numerical simulations (Goodwin & Whitworth 2004), although some coherence in the 
initial velocity dispersion is required. Another explanation is that this difference is a 
consequence of a more clustered distribution of the densest gas from which stars form at 
the smallest spatial scales in the molecular cloud complexes, according to a multifractal 
scenario (Chappell & Scalo 2001). Maybe the star formation process itself modifies in 
some (unknown) way the underlying geometry generating distributions of stars that can 
be very different from the distribution of gas in the parental clouds. Finally, one pos- 
sibility is that the fractal dimension of the interstellar medium in the Galaxy does not 
have a universal value and therefore some regions form stars distributed following more 
clustered patterns. There is no a priori reason for assuming that Df has nearly the same 
value everywhere in the Galaxy independently of either the dominant physical processes 
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Figure 2. Structure parameter Q as a func- 
tion of the logarithm of age divided by the 
tidal radius, which is nearly proportional to 
age in crossing times units. The dashed line at 
Q = 0.8 roughly separates radial from fractal 
clustering. The best linear fit (equation at the 
top) is represented by a solid line. 



Figure 3. Calculated correlation dimension 
as a function of age (in crossing time units). 
The best linear fit (equation at the top) is 
represented by a solid line. As reference, hor- 
izontal dashed lines indicate the values cor- 
responding to three-dimensional distributions 
with fractal dimensions of Df — 2.0 and 2.5. 



or environmental variables. Recent simulations of supersonic isothermal turbulence done 
by Federrath et al. (2009) showed that compressive forcing yields fractal dimension val- 
ues for the interstellar medium significantly smaller (Df ~ 2.3) compared to solenoidal 
forcing (Df ~ 2.6). Thus, Df could be very different from region to region in the Galaxy 
depending on the main physical processes driving the turbulence. At least at galactic 
scales, it has been shown that there are significant differences in the fractal dimension 
of the distribution of star forming sites among the galaxies, contrary to the universal 
picture previously claimed in the literature (see Sanchez & Alfaro 2008). So that the 
possibility of a non-universal fractal dimension for the interstellar medium in the Galaxy 
cannot, in principle, be ruled out. 
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